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Abstract. We discuss the theoretical underpinnings for the extraordinary success of perturbative 
QCD in the description of HERA data. In particular, we examine recent progress in the understand- 
ing of perturbative QCD at small x. We explain the relation between evolution equations in Q 2 and 
x, and how they can be used for simultaneous resummation of the relevant large logs at HERA. We 
show that while the HERA data can be understood within our current knowledge of the perturbative 
expansion of QCD, they pose stringent constraints on the perturbatively inaccessible behaviour of 
QCD in the Regge limit. 
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PERTURBATIVE QCD AT HERA 

QCD has been tested at HERA [1,2] over the last several years to an accuracy which 
is now comparable to that of tests of the electroweak sector at LEP: perturbative QCD 
turns out to provide an embarrassingly successful description of the HERA data, even 
in kinematic regions where simple fixed-order perturbative predictions should fail. This 
success is most strikingly demonstrated by the comparison with the data of the scaling 
violations of structure functions predicted by the QCD evolution equations [3, 4]: the 
data agree with the theory over five orders of magnitude in both x and Q 2 . 

The significance of this sort of result is somewhat obscured by the need to fit the 
shape of parton distributions at a reference scale, which might suggest that deviations 
from the predicted behaviour could be accommodated by changing the shape of the 
parton distribution. However, this is not true because of the predictive nature of the 
QCD result: given the shape of partons at one scale, there is no freedom left to fit the 
data at other scales. This predictivity is particularly transparent in the small x region, 
where the fixed-order QCD result actually becomes asymptotically independent of the 
parton distribution, apart from an overall normalization. Indeed, the data for lni<2 plotted 

versus the variable a = In -f In - are predicted to lie on a straight line, with universal 

slope 2y = 12/ ^33 — 2m} (double asymptotic scaling [5, 6]). The predicted scaling is 
spectacularly borne out by the data, as shown in fig. 1: in fact, the data are now so 
accurate that one can see the change in slope when passing the b threshold, and indeed 
double scaling is only manifest if one separates data in the regions where a y runs with 
Nf = 4 from those with Nf = 5. 2 Equally good agreement with fixed-order perturbation 
theory is seen when considering less inclusive observables. 

This agreement of the data with fixed-order perturbative QCD computations is very 
surprising, in that the perturbative expansion receives contributions of order CXyln^ so 
one would expect higher-order corrections to be non-negligible whenever oc v In i ^ 1 , 
i.e. in most of the HERA region. As is well known, the resummation of leading In i 
(LLx) contributions to gluon-gluon scattering, and thus to a wide class of hard pro- 
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2 The fact that the observed slope is somewhat smaller than the predicted one, especially at low Q 2 , is due 
to NLO corrections [7] as well as corrections due to the "small" eigenvalue of perturbative evolution [8]. 
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FIGURE 1. Double asymptotic scaling of the HI data [4]. The scaling variable o = 
ln(xo/x)ln(a y (go)/c.v(6 2 )) is defined with xq =0.1, go = 1 GeV; therescaling factor/?/? is as in Ref. [5]. 
Only data with p > 1, o > 1, Q 2 > 4 GeV 2 ; x < 0.03 are plotted. Left: Q 2 < m 2 ; right: Q 2 > m\. The 
straight line is the asymptotic prediction. 



cesses, including small x scaling violations of structure functions, is accomplished by 
means of the BFKL evolution equation [9, 10, 11]. Matching the BFKL approach to 
standard perturbative computation, however, is nontrivial [12, 13], while the BFKL 
equation itself seems to be unstable towards the inclusion of higher order correc- 
tions [14]. Hence, the main problem in understanding HERA physics, i.e. perturbative 
QCD at small x is that of establishing "consistency of the BFKL approach with the 
more standard DGLAP [15, 16] evolution equations" [17], which embody the leading 
In Q 2 (LLQ 2 ) resummation on which perturbative QCD is based. This problem is now 
solved [18, 19, 20], and on the basis of this solution it is possible to combine the avail- 
able information on perturbation theory at small x, and use it to explain the unexpected 
success of fixed-order calculations. 



Let us for definiteness consider the prototype problem of the description of small x 
scaling violations of parton distributions. For simplicity, consider the case of a single 
parton distribution G(x, Q 2 ), which can be thought of as the dominant eigenvector 
of perturbative evolution. Scaling violations are then described by the Altarelli-Parisi 
equation satisfied by G{x, Q 2 ), and thus summarized by the Altarelli-Parisi splitting 
function P(x, a s ) [15]. 

The basic result which allows the study of scaling violations at small x is duality 
of perturbative evolution [21, 19, 22], namely, the fact that, because the Altarelli- 
Parisi equation is an integro-differential equation in the two variables t = ln<2 2 /A 2 and 
£ = l/x, it can be equivalently cast in the form of a differential equation in t satisfied by 
the jc-Mellin transform 



DUALITY 




(1) 



or a differential equation in £ satisfied by the <2 2 -Mellin transform 
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of the parton distribution. The pair of dual evolution equations are 

jG(N,t) = y(N,a s )G(N,t) (3) 
^G&M) = x(M,a s )G(^,M), (4) 

where eq. (3) is the standard renormalization-group equation, with anomalous dimen- 
sion y(N,t), and eq. (4) is essentially the BFKL equation. Duality is the statement that 
the solutions of these two equations coincide to all perturbative orders, up to power 
suppressed corrections, provided their kernels are related by 

X(y(N,a s ),a s )=N. (5) 

This means that the BFKL and Altarelli-Parisi equations describe the same physics: 
it is the choice of the kernel to be used in the evolution equation which determines 
which is the large scale which is resummed. We can then discuss the construction and 
resummation of the kernel irrespective of the specific evolution equation where it is 
used, with the understanding that the kernel can be equivalently viewed as a y(N, a s ) or 
a %(M,a s ), the two being related by eq. (5). Before doing this, we sketch how duality 
can be proven order by order in perturbation theory. 



Fixed coupling 

Perturbative duality is most easy to prove when the coupling does not run, since in 
this case the two scales t and £, appear in the Altarelli-Parisi equation in a completely 
symmetric way. It is convenient to introduce the double-Mellin transform G{N,M) of 
the parton distribution. The solution to the Altarelli-Parisi equation in M,N space has 
the form (which can be e.g. obtained by performing an M-mellin transform eq. (2) of 
the solution to the renormalization-group eq. (3)) 

where Gq(N) is a boundary condition at a reference scale /a 2 . 

The inverse Mellin transform of eq. (6) coincides with the residue of the simple pole 
in the M plane of e tM G(N,M), and thus its scale dependence is entirely determined by 
the location of the simple pole of G(N 7 M) (6) , namely, the solution to the equation 



M = y(N,a s ). 



(7) 



The pole condition Eq. (7) can be equivalently viewed as an implicit equation for N: 
N = x(M, cCy), where % is related to yby eq. (5). Hence, the function 



corresponds to the same G{t,x) as eq.(6), because the location of the respective poles in 
the M plane are the same, while the residues are also the same, provided the boundary 
conditions are matched by 

Go(A ° = -fSM))- (9) 

Eq. (8) is immediately recognized as the /V-Mellin of the solution to the evolution 
equation (8) with boundary condition Fq{M) (at some reference x = xo), which is what 
we set out to prove. In general, the analytic continuation of the function % defined 
by eq. (5) will be such that eq. (7) has more than one solution (i.e. y is multivalued). 
In this case, poles further to the left in the M plane correspond to power-suppressed 
contributions, while poles to the right correspond to contributions beyond perturbation 
theory (they do not contribute when the inverse M-Mellin integral is computed along 
the integration path which corresponds to the perturbative region). 

It is easy to see that upon duality the leading-order % = CC v %o is mapped onto the 
leading singular y= y s (a s /N), and conversely the leading-order y= a s yo is mapped 
onto the leading singular % = % s (a s /M). In general, the expansion of % in powers of 
a s at fixed M is mapped onto the expansion of y in powers of cc^ at fixed OL s /N, and 
conversely. So in particular at LLQ 2 it is enough to consider yo or % s , and at LLx it is 
enough to consider y s or %o- The running of the coupling is a LLQ 2 but NLLx effect, so 
beyond LLx the discussion given so far is insufficient. 



Running coupling 

The generalization of duality to the running coupling case is nontrivial because the 
running of the coupling breaks the symmetry of the two scales £, and t in the Altarelli- 
Parisi equation. Indeed, upon M-Mellin transform (2) the usual one-loop running cou- 
pling becomes the differential operator 

where da s /dt = — PoOC 2 . 

Consider for simplicity the LLx x-evolution equation, i.e. eq. (4) with % = a s %o(M), 
and include running coupling effects by replacing a s with the differential operator 
eq. (10). We can solve the equation perturbatively by expanding the solution in powers 
of at fixed CL S /N: the leading-order solution is given by eq. (8), the next-to-leading 
order is obtained by substituting this back into the equation and retaining terms up to 



order Potty, and so on [18]. We can then determine the associate G(N,t) by inverting 
the M-Mellin, and try to see whether this G(N, t) could be obtained as the solution of a 
renormalization group (RG) equation (3). 

The inverse Mellin is again given by the residue of the pole of e tM G(N 7 M) in the 
M-plane, where G(N,M) is now the perturbative solution. When trying to identify this 
with a solution to eq. (3) there are two potential sources of trouble. The first is that now 
the perturbative solution at order (a s $o) n has a (2n+ l)-st order pole. Therefore, the 
scale-dependence of the inverse Mellin is now a function of both oc s and t, whereas 
the solution of a RG equation depends on t only through the running of oc v . Hence it is 
not obvious that a dual anomalous dimension will exist at all. The second is that even 
if a dual y does exist, it is not obvious that it will depend only on % and not also on 
the boundary condition Fq(M) eq. (8): in such case, the running of the coupling in the 
^-evolution equation would entail a breaking of factorization. 

However, explicit calculation shows that it is possible to match the anomalous dimen- 
sion and the boundary condition order by order in perturbation theory in such a way that 
both duality and factorization are respected. Namely, the solution to the leading-twist 
running coupling ^-evolution eq. (4) with kernel a s %o and boundary condition Gq(M) 
is the same as that of the renormalization group eq. (3) with boundary conditions and 
anomalous dimension given by 

y(a s (t),a s (t)/N) = y s (a s (t)/N) + a s (t)$ Ay ss (a s (t)/N) + 

+(a. s .(0(3o) 2 Ay,,v(a.,(0/^V) +0(a.,(OPo) 3 (H) 
G (a,,AO = G (N) +a$ A^G (N) + (a,p ) 2 A (2) G (^V) +O(a,p ) 3 , (12) 

where the leading terms y s and Go(N) are given by eqs. (5) and (9) respectively. The 
subleading corrections are 

AY.,- = -|| (13) 

A (i) Go(iV) = ^Po-'XfiW-M') (14) 

where all derivatives are with respect to the arguments of %o(M) and Fq(M), which 
are then evaluated as functions of y s (a s /N). The sub-subleading correction to the 
anomalous dimension is 

1 W /3 _ 1 f.y> y"y"' I 2y' 2 y"" 

a». _ „2 13 >Co 10 XqXqXo + ^Iq Xp ni -s 

KYsss — — Xo ZTTs ' v 1 3 ) 

24 Xo 

and we omit the very lengthy expression for A^Gq(N). The fact that duality and 
factorization hold up to NNLLx is nontrivial, and suggests that they should hold to all 
orders. An all-order proof can be in fact constructed [23]. 

Once the corrections to duality eq. (12) are determined, they can be formally re- 
interpreted as additional contributions to %: namely, one can impose that the duality 
eq. (5) be respected, in which case the kernel to be used in it is an "effective" %, obtained 



from the kernel of the ^-evolution eq. (4) by adding to it running coupling corrections 
order by order in perturbation theory: %o will be free of such correction, %i will receive 
a correction 

A%i = p ^ -j— — , (16) 

and so forth. Applying duality to the known one-loop anomalous dimensions yo thus 
gives us the resummation of the all-order singular contributions %(a s /M) to this ef- 
fective %, which include the running coupling correction eq. (12) and its higher-order 
generalizations. 



RESUMMATION 

Because the first two orders of the expansion of % in powers of a s at fixed M and of 
the expansion of y in powers of cc s at fixed ./V are known, it is possible to exploit duality 
of perturbative evolution to combine this information into anomalous dimension which 
accomplish the simultaneous resummation of leading and next-to-leading logs of x and 
Q 2 . In fact, it turns out that both a small M and a small N resummation of anomalous 
dimensions are necessary in order to obtain a stable perturbative expansion, while 
unresummed anomalous dimensions leads to instabilities. Both sources of instability are 
generic consequences of the structure of the perturbative expansion, and could have been 
predicted before the actual explicit computation [14] of subleading small-* corrections. 



Small M 

The perturbative expansion of % at fixed M is very badly behaved in the vicinity of 
M ~ 0: at M = 0, %o nas a simple pole, %i has a double pole and so on. In practice, this 
spoils the behaviour of % in most of the physical region < M < 1. Because l/M k is 

the Mellin transform of ln^ 1 (Q 2 /A 2 ), these singularities correspond to logs of Q 2 
which are left unresummed in a LLx or NLLx approach [24]. 

The resummation of these contributions may be understood in terms of momentum 
conservation, which implies that y(l,a. v ) = (note our definition of the /V-Mellin 
transform (1), and also that y is to be identified with the large eigenvector of the 
anomalous dimension matrix). The duality eq. (5) then implies that a momentum- 
conserving % must satisfy %(0, a s ) = 1. This, together with the requirement that % admits 
a perturbative expansion in powers of a s , implies that in the vicinity of M = 0, the 
generic behaviour of the kernel is 

~ - - °^ I °^ I (17) 

m^ooc, + kM kM (kM) 2 (kM) 3 "' ' 

where K is a numerical constant which turns out to be K = tz/Ca- Hence we understand 
that there must be an alternating-sign series of poles at M — 0, which sums up to a 
regular behaviour. In fact, we can systematically resum singular contributions to % to 



all orders in oc s by including in % the terms % s (a s /M) derived from the leading order 
Jo(N), and similarly at next-to-leading order, and so on. Because the usual anomalous 
dimension automatically respects momentum conservation order by order in a s , in order 
to remove the small M instability of the expansion of % at fixed M, it is sufficient to 
improve the expansion by promoting it to a "double leading" expansion which combines 
the expansions in powers of 0C S at fixed M and at fixed CL S /M [19]. For example, 
at leading order % = a s %o(M) + % s {cl s /M) — d.c, where the subtraction refers to the 
double-counting of the a s /M term which is present both in a v %o and in % S {(X S /M) . This 
expansion of % is dual eq. (5) to an analogous expansion of J, where at leading order 
Y= OCvYo(M) +y s (a s /M) —d.c, and so forth. Both expansions are well behaved at small 
M, i.e. large N. At this level, it is already clear that the impact of the inclusion of small-* 
corrections is moderate: indeed, it turns out that the double-leading kernel is quite close 
to the usual two-loop kernel, except at the smallest values of N, i.e. in the neighbourhood 
of the minimum of %(M) [19]. 



Small N 

The improved double-leading expansion of the anomalous dimension still requires 
resummation at small N. This is because, even though the next-to-leading correction 
to the double-leading evolution kernel is small for all fixed M, it is actually large if 
N is fixed and small. This in turn follow from the fact that the leading % kernel has a 
minimum, so the small N = % region corresponds by duality eq. (5) to the vicinity of the 
minimum where the kernel is almost parallel to the y = M axis. 

At small N, unlike at small M, there is no principle like momentum conservation 
which may provide a fixed point of the expansion and thus fix the all-order behaviour. 
The only way out is thus to treat this all-order behaviour as a free parameter. Namely, 
we introduce a parameter X which is equal to the value of the all-order kernel % at its 
minimum, and then we expand about this all-order minimum. In practice, this means 
that we reorganize the expansion of % according to [18] 

%{M,a s ) = a s % (M) + a 2 s Xi{M) + ... 

= a s % (M) + a 2 MM) + ..., (18) 

where 

oo 

a s %o(M, a,) = a,%o(M) + £ a? +1 c„, fc(M) = Xi(M) - c u (19) 

n=l 

and the constants c,- are chosen in such a way that 

X = a s jfo(£) = cWCo(2)+AX. ( 2 °) 

is the all-order minimum of %. Of course, in practice phenomenological predictions will 
only be sensitive to the value of X in the region where very small values of Af are probed, 
i.e. at very small x. 
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FIGURE 2. x 2 (left) and starting gluon slope G(x,4 GeV 2 ) - (right) for the fit [22] to the 95 HI 
data [4] as a function of the resummation parameter X eq. (20), for the two resummation prescriptions 
discussed in text. The fits are performed with tt s (M z ) = 0.1 19. 



PHENOMENOLOGY 

Using duality and the resummation discussed above, one can construct resummed ex- 
pressions for anomalous dimensions and coefficient functions, and wind up with re- 
summed expressions for physical observables which may be directly compared to the 
data. The need to resum the small N behaviour entails that phenomenological predictions 
will necessarily depend on the parameter X eq. (20). When the resummed double-leading 
expansion is constructed, a further ambiguity arises in the treatment of double-counting 
terms. This ambiguity is related to the nature of the small N singularities of the anoma- 
lous dimension, which control the asymptotic small x behaviour. Specifically, according 
to the way the double-counting is treated, the N = poles of the one- and two-loop re- 
sult may survive in the resummed result ('S-resummation') or not ('R-resummation'). 
Both alternatives are compatible with the known low-order information on the evolution 
kernel, and can be taken as two extreme resummation schemes which parametrize our 
ignorance of higher order perturbative terms. Since the resummed terms also have a cut 
starting at N = X, whether or not these low-Af poles are present only makes a difference 
if X turns out to be small, X & 0.3. 

The % 2 and starting gluon slope for a fit [22] to the recent HI data [4] for the deep- 
inelastic cross section are shown in figure 2, as a function of X and for the two different 
resummation prescriptions. It is clear that if the perturbative N = poles do not survive 
the resummation (R resummation) then only a fine-tuned value of X 0.2 is acceptable, 
whereas if they do survive (S resummation) essentially any X ^ gives a good fit. 

Figure 2 demonstrates that it is possible to accommodate the success of simple fixed- 
order approach within a fully resummed scheme, and in fact the resummed calculation 
is in somewhat better agreement with the data than the fixed order one. Even though the 
effects of the resummation are necessarily small (otherwise the success of the fixed order 
prediction could not be explained) they do have a significant impact in the extraction of 



the parton distribution: the gluon comes out to be significantly more valence-like than 
in an unresummed fit. Hence, the use of resummed perturbation theory is crucial for the 
extraction of reliable parton distributions at small x. 

From a theoretical point of view, we see that current data already pose very stringent 
constraints on the unknown high-orders of the perturbative expansion: only a rather 
soft high-energy behaviour of the deep-inelastic cross-section is compatible with the 
data. Further progress in the understanding of the Regge limit is likely to require 
either genuinely nonperturbative input, or an extension of the standard perturbative 
domain [18]. 
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